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Abstract. Based on a recent classification of subalgebras of the symmetry 
algebra of the Zabolotskaya-Khokhlov equation, all similarity reductions of 
this equation into ordinary differential equations are obtained. Large classes 
of group invariant solutions of the equation are also determined, and some 
properties of these solutions are discussed. 



1. Introduction 

The Zabolotskaya-Khokhlov equation, which is a sound wave propagation model 
derived from the incompressible Navier-Stokes equation [TJ [2J [5] has been studied 
from the Lie group approach in a number of papers in which similarity [3 El [7] and 
other types of reductions [10] were obtained. A variant of different order of this 
equation was also discussed in [4]. However, as far as reductions to ordinary differ- 
ential equations (ODEs) are concerned, all the results obtained in these papers were 
severely limited by a number of anzatz, and generally by the lack of classification 
into subalgebras of the symmetry algebra of the equation. Although a complete 
group classification of the symmetry algebra for the (2 + l)-dimcnsional version of 
this equation was recently given in [7j, this only led, due to space limitations, to 
the determination of all similarity reductions to (1 + l)-dimensional models of the 
equation in that paper. 

In this paper, we consider the (2 + l)-dimensional version of the Zabolotskaya- 
Khokhlov (ZK) equation, which is given by 

A(t, x, y, u) = u xt - {uu x ) x - u yy = 0. (1.1) 

Using the complete classification of two-dimensional subalgebras of the symmetry 
algebra of this equation recently proposed in [7], all similarity reductions of this 
equation to ODEs are obtained. These ODEs have the advantage of being much 
easier to solve compared the (1 -I- l)-dimensional models obtained in similarity 
reductions by one-dimensional subalgebras, and thus new exact similarity solutions 
are obtained, and some of their properties are discussed. It is shown in particular 
that all nonlinear reduced ODEs are non linearizable. 



2. Classification of two-dimensional subalgebras 

We wish to recall in this section a result on the classification of two-dimensional 
subalgebras of the symmetry algebra of the ZK equation obtained in a recent pa- 
per [7] . The Lie symmetry algebra L of the ZK equation itself is well-known [9] . 
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It is an infinite dimensional algebra with generators 

v = 2xd x + yd y + 2ud u (2.1a) 
x(s) = gd x - g'du (2.1b) 

y(h) = lyh'd x + hd y -^yh"d u (2.1c) 

<f) = fd t + \ (2xf + y 2 f") d x + ^f'd v + i (-4uf 2xf" - y 2 f">) d u 

(2.1d) 

where /, g, h are arbitrary functions of the time variable t 1 and where a prime 
represents a derivative with respect to t. It is easy to see that for the Lie algebra 
L, the commutation relations are given by 

[v ,x( 3 )] = -2x( ff ), [x( S ),y(/i)] = (2.2a) 

[v ,y(M] = -y(h), [x(j),«(/)] - x(/V3-/.g') (2.2b) 

[v ,z(/)] = o, [y(h)Xf)} = y(fA-A') (2.2c) 

[x( 5l ),x(.g 2 )] = [y(/*i),y(M] - x((h 1 ti 2 -h' 1 h 2 )/2) (2.2d) 

[z(/i),z(/ 2 )] = z(/ x ^-/(/ 2 ). (2.2e) 

These commutation relations show that a Kac-Moody-Virasoro (kmv) structure 
can be associated with an infinite dimensional subalgebra of L [7], and the latter 
property tends to associate with integrability the (2 + l)-dimensional ZK equa- 
tion [T^l E3 Q3] , which, as is well-known, can be linearized by a generalized hodo- 
graph transformation [16] . 

In order to classify subalgebras of L under the adjoint action of its Lie group 
Q, we need to have an explicit expression for w(e) = Ad(exp(ev))wo, for every 
pair of generators v, Wo of L. However, using the commutation relations (|2.2p such 
an expression can easily be obtained either by interpreting w(e) as the flow of Ad 
through wo of the one-parameter subgroup generated by v, or again by rewriting 
w(e) in terms of the Lie series (see [TTJ P. 205]). The required classification of sub- 
algebras of L under the adjoint action of Q can henceforth be achieved by applying 
known techniques |15l 111] . In this way, all one-dimensional and two-dimensional 
subalgebras of L were classified in a recent paper |7] and the corresponding list of 
canonical forms of non-equivalent two-dimensional subalgebras of L can be repre- 
sented in terms of their generators as follows, where fcrj, k\, c\, c%, and C3 are free 
parameters, unless otherwise stated. 
(1) Abelian subalgebras 

-^2,1 = {v , fcovo + z (1) } 

^2,2 = {A:iv + z (1) , /c v + X( Cie 2fc lt) + y (c2£ fc lt ) + z (c3 )} 
^2,3 = {x ( i 3 ,x ff }, ( S V0) 
^2,4 = {x (1) ,y h } 

^2,5 = {x(i) , fc v + z (6fcot+C3) } , k 2 Q + c\ ^ 

= {y(i),x g +y (c2) } 
^2,7 = (y(i) , fcovo + y (c2 ) + z (( _ 3/ 2)( feot+C 3)) } , k\ + cj ^ 
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(2) Non — abelian subalgebras 

■#2,8 = {z (1) ,fc v +X( Cl ) +Y( C2 ) +Z( t+C3 )} 

-%,9 = {x (1 ),V } 

-^2,10 = {x (1) ,fc v +z ((3 _ 6A , o)t+C3) } , (fc -5) 2 + c^O 

-22,ii = {y(i),v +y (c2) } 

-22,12 = {y(l) , fcoVo + y( C2 ) + Z((3/2)(l-fc )t+c 3 ) } > ( fc ~ I) 2 + C§ ^ 

In this list, «5f 2 ,j represents the jth canonical form of a two-dimensional subal- 
gebra of L. 

3. Similarity reductions to ODEs 

Similarity reductions of the ZK equation by symmetry subgroups whose actions 
are semi-regular with orbits of dimension s, where < s < p, and where p = 3 is 
the number of independent variables, will yield equations in s fewer independent 
variables. Each of the canonical forms J>?2,j of two dimensional subalgebras of 
L turns out to have two dimensional orbits, and hence allow a reduction of the 
equation into an ODE, provided certain regularity conditions are satisfied. The 
usual procedure for this reduction is well known and consists in this case in finding 
a pair of invariants z — z(t,x,y,u) and w — w(t,x,y,u) of the corresponding 
subgroup action and to consider one of them, say z, as the new independent variable 
and the other one as the new dependent variable for the reduced equation. 

To find the required similarity reductions, we shall treat each of the twelve 
canonical forms _£? 2 ,j separately. In each case, we will indicate a reduction formula 
which will consist in an expression for the solution u of the ZK equation in terms 
of w = w(z), and an expression for z, and then the corresponding reduced equa- 
tions. The parameters ko, fci, Ci, C2, C3 are the same as those appearing in the list 
of classification of canonical forms -22,j given in the previous section. 

3.1. Reduction by abelian subalgebras. 

3.1.1. Reduction by _£?2,i- The reduction formula is 

u~xw(z), z~y 2 /x, (3-1) 

and the reduced equation is 

- w 2 + 2{wz - l)io' - z 2 w' 2 - (4z + wz 2 )w" = 0. (3.2) 

3.1.2. Reduction by J?2,2- The reduction formula is 

(4c lC2 e tfcl y + 2c iy 2 (k - c 3 fci) + cK-4e + y 2 ^)) 

2(c 2 e tfc i +y{k -c 3 k 1 )) 1 ' &) 

-yk x {c\yk\ + 4ci(2c 2 e tfcl + y{k - c 3 ki))) + 4(c 2 e tfcl + y{k - c 3 fci)) 2 w(z) 

4c 2 

(3.3b) 

and the reduced equation which includes the case c 2 — is given by 
- 8wz 3 (fc - c 3 fci) 2 + 2z 3 k 1 (4:C 1 koclk 1 - 4cic 3 fci) 

+ (c 2 + c\w + 8z 4 {k - c 3 /ci) 2 + 3z 2 (c^i + 2ci(fc - c 3 fci)))«;' - c\zw 2 (3.4) 
+ ({-c\ - c\w)z - 4z 5 (fc - c 3 fci) 2 - 2z 3 (2 Cl fc + c\k x - 2c 1 c z k 1 )) w" = 
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3.1.3. Reduction by _£?2,3- In this case, the similarity variables z and w are given by 
z = t, w = y. Hence the corresponding Lie point transformation does not satisfy 
the transversality condition, since the equation w(z) = y cannot be solved for u, 
and there is no reduced equation in this case. 

3.1.4. Reduction by J^2,i- The reduction formula is 

z = t, u = (w(z) - y 2 )h"/Ah. (3.5) 

In this case the reduced equation degenerates simply into the condition h" = 0, 
which in turn yields only the trivial solution u = 0. 

3.1.5. Reduction by J^2,b- The reduction formula is 

(c 3 + 6M) 5 w{z) . . 

z = —t — ' u= (c 3 + 6k ty^y forfc( ^° (3 ' 6a) 

z = y, u = w(z), for fc = (3.6b) 
and the corresponding reduced equation is 

7w' + 6zw" = 0, for fc ^ (3.7a) 

w" = 0, for fc = 0. (3.7b) 

3.1.6. Reduction by ££<i$. The reduction formula is 

z = t, u = g'(w(z) - x/g). (3.8) 

Here again the reduced equation is simply a degenerated condition of the form 
g" = 0, which gives rise to a solution of the form 



„ = -_j (3.9) 

where a and (3 are arbitrary constants, while w(t) is an arbitrary function of t. 

3.1.7. Reduction by _£?2,7- The reduction formula is 

z = x(c 3 + k t), u = w(z)/(c 3 + k t) 2 , + (3.10) 

and this leads to the reduced equation 

k Q w' + w' 2 + (w - zk a )w" = 0. (3.11) 

3.2. Reduction by non abelian subalgebras. 

3.2.1. Reduction by _£?2,8- The reduction formula in this case is 

(6fc + l)x + 3ci ,,„, „. „ .. , . /, . 1 2 

u = ((6feo + l)x + 3 Cl ) r w(z), fc ^ 



((3fc + 2) ? y + 3c 2 ) 1 ' " / v n y » / 6 . 3 

^=^f£F' (fco * -1/6) (3.12a) 

z = (ci -x)e y l c \ u= { Cl - x) 2 w{z), (fe = -2/3) (3.12b) 

where s = (6fc + l)/(3fc + 2) and r = (6k — 2)/(6fc + l). This leads to the following 
reduced equation when ko ^ —1/6 and kg ^ —2/3. 

- r(-l + 2r)w 2 z r {\ + 6fc ) 2 + (-«(1 + s)z 3 (2 + 3fc ) 2 - 4rwz r+1 (l + 6k Q ) 2 )w' 

- z r+2 (l + 6fc ) V 2 + (s 2 z 4 (2 + 3fc ) 2 - wz r+2 (l + 6k ) 2 )w" = 0. 

(3.13) 



The other reduced equations for fco = —1/6 and fco = —2/3 are respectively given 

by 

-2w 2 + (3wz - §c\z 2 )w' - z 2 w' 2 + {~wz 2 - Ac\z z )w" = 0, (3.14a) 
-6clw 2 -{z + 8c 2 2 wz)w' - c 2 2 z 2 w 12 - (z 2 + c\wz 2 )w" = 0, (3.14b) 

3.2.2. Reduction by Jz?2,9- The reduction formula is 

z = t, u = y 2 w(z), (3.15) 
which only results in the trivial equation w = 0, and solution u = 0. 

3.2.3. Reduction by J^2,w- The reduction formula is 

z= C3 + (3 y 7 u=(c 3 + (3-6k )t) s w{z), for fc ^1/2 (3.16a) 

z = e (-t/2cs) yj „ = e t/ca w ( z ^ for fc = 1/2. (3.16b) 

where r = (3 — 6fco)/(2 — 3fco) an d s = (6fco — 2)/(3 — 6fco)7 an d the corresponding 
reduced equation is 

(-5 + 9fco)u/ + 3z(-l + 2fco)w" = 0, for fc 7^ 1/2 (3.17a) 

w" = 0, for fc = 1/2 (3.17b) 

3.2.4. Reduction by JSf2.11- The reduction formula is 

2 = i, u = 2:10(2), (3.18) 
and this gives rise to the reduced equation 

w 2 ~w' = Q. (3.19) 

3.2.5. Reduction by Jzf2,i2- The reduction formula is 

z = (3(l-fc )t + 2c 3 )x s , u= (3(1 - fc )t + 2c 3 ) r u;(2), for fc ^ 1,-1/3 

2 = e (2 * /c3) /x, u = e (2t/c3) w(z), for fc = 1 

z = x, u = w(z)/(c3 + 2i), for fc = — 1/3 

where s = 3(fc - 1)/(1 + 3fc ) and r = (2(3fc - 1))/(3(1 - fc )). This yields for 
fc ^ 1,-1/3, the reduced equation 

((-1 - s)wz 1+r + 3(1 + r)(fc - 1)) w' - sz 2+r w' 2 

+ (swz 2+r + 3z(fc - l))w" = 0, (3 ' 21) 

while the other reduced equations are 

(4 + 2c 3 wz)w' + c 3 z 2 w' 2 + z(2 + c 3 wz)w" = 0, for fc = 1 (3.22a) 

2w' + w' 2 + ww" = 0, for fc = -1/3 (3.22b) 

All of the reduced ODEs that we have thus obtained for the ZK equation belong 
to a more general family of differential equations of the form 

Ai +A 2 w' + A 3 w' 2 + A A w" = (3.23) 
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where the Aj are polynomials in w whose coefficients are functions of z, and of the 
form 

Ai = aiw 2 + a 2 wz r + a 3 wz 3 + a^z 3 (3.24a) 
A 2 = (3 lW z r+1 + fowz + (3 3 w + P 4 (z) (3.24b) 
A 3 = 7l z r+2 + P 2 (z) (3.24c) 
A A = 5 1 wz r+2 + 6 2 wz 2 + S 3 wz + S 4 w + P 5 (z), (3.24d) 

where the a j , /3j , 7j , <5j and r are free constants, while P n {z) represents a polynomial 
of degree n in z. These equations are generally nonlinear, except in the case of 
reductions by J2? 2 ,5 and -S?2,io where the corresponding equations of the form (|3.23|) 
turn out to be genuine linear equations of the simple form 

aiw' + a 2 zw" = 0, (3.25) 

where and a 2 are constants. 



4. Exact solutions and their properties 

Although explicit expressions for solutions of the ZK Equation (jl.ip can be ob- 
tained from the reduced equations of the previous section, we shall also be interested 
in the physical relevance of these solutions. In the ZK model p], physically rele- 
vant solutions must be bounded functions of the spatial variable x, aligned with 
the primary direction of propagation, and the time-like variable t. 

In this sense, not all solutions of the ZK equation are physically relevant, and 
this is the case with the solution given for example by ()3.9() . It is also easy to verify 
that for reduced ODEs which are linear, relevant solutions occur only in the case of 
Eq. ()3.7|) . while none occurs for the other general case given by (|3.17|) . Thanks to 
the reduction formulas, and the fact that the integration of linear equations of the 
form (|3.25p poses no problem, explicit expressions of solutions are not necessary in 
the case of linear reduced equations. 

For nonlinear reduced equations, the simplest and most frequent equations of 
the form (|3.23j) are equations in which the coefficients Aj are simple polynomial 
functions in w and z. More precisely, these are given by Ai = aw 2 , and A 3 = (3z 2 , 
where a and (3 are constants, while A 2 and A4 are linear in w. We call Type A this 
class of simpler reduced nonlinear ODEs, and Type B its complement among the 
nonlinear reduced ODEs. Type A equations often correspond to subcases of Type 
B equations which in turn correspond to general cases of reductions by -22,2,-^2,8 
and JS?2,i2- We can partition Type A into Type A\ consisting of equations admitting 
only one symmetry, and Type A 2 , consisting of equations admitting two or more 
symmetries. All Type A equations can generally be reduced to an Abel equation 
of the first kind, which as is well-known are generally difficult to solve [17l ITS] . 
However, when any such equation has two known symmetries, it is possible to 
obtain solutions of the nonlinear equation at least in an implicit or parametric 
form. It will be sufficient to illustrate this point by treating a number of examples. 

If we consider for instance Type A reduced equation (|3 - 2[) corresponding to 
the subalgebra -£?2.i, we observe that this equation has only a one-dimensional 
symmetry algebra generated by 



v = z d z — w <9, 



(4.1) 
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with rectifying coordinates r = wz and W = — ln(w). In terms of these new vari- 
ables, Eq. p.2p takes the form 

- 1 + (2 - 5r)W' + (-8r + \\r 2 )W' 2 - 6r 2 (l + r)W' 3 + r(4 + r)iy" = 0, (4.2) 

which is clearly an equation of Abel of the first kind for Z — W' . However, even 
the transformed equation (|4.2|) also has only one symmetry, and thus it is difficult 
to find its solutions by Lie groups methods. Similar results apply to all other Type 
Ai reduced ODEs such as (|3"TTj) and ([3~T4"j) . 

Among Type A2 equations the simplest is Eq. (|3.19p . which is of the first order 
and corresponds to a reduction by -S?2,ii> but the corresponding solution is of no 
physical interest as already mentioned. In all other Type A 2 cases, the equation 
is of order two and can also be reduced to an Abel equation of the first kind. A 
simple example for this case is Eq. (|3.22b[l . It admits the symmetry v = d z , in 
terms of whose rectifying variables r = w,W = z, the equation takes the form 

2W' 2 + W'-rW" = 0, (4.3) 

and thus reduces to a Riccati equation for Z = W', which is just a degenerated 
form of an Abel equation. However, Eq. (|3.22b[) can be solved directly, and gives 
rise to a a bounded solution 

A 

u = —, r 

2(c 3 + 2t) 

where A and B are arbitrary constants, and where ProductLog(z) is the principal 
solution for w in the equation z = we w . Equation (|3.22a[) has the same properties, 
although it is less obvious to solve. It has two symmetries v = z 2 d z + (2/C3) d w and 
w = z d z — w d w satisfy a commutation relation of the form [w, v] = v. Thus by 
the well- know procedure of [llj , the equation should first be reduced by the normal 
symmetry subalgebra generated by v. Before trying to solve (|3.22a|l . we notice that 
in terms of the rectifying coordinates r — 2/z + c^w and W — 3csw/2 + 2/z of v, 
Eq. (I3.22aj) takes the form 

--— W' + — W' 2 --W' 3 + W" = 0, (4.5) 
r r r r 

which is clearly an Abel equation for Z — W' . On the other hand, in terms of the 
invariants £ = 2/z + C3W and X — z 2 w' of the first prolongation of v, Eq. (j3.22a[) 
reduces to the first order equation 

- c 3 X 2 + f(2 - c 3 X)X' = 0. (4.6) 

In terms of the new variables £ and X, w reduces in a quotient space to a symme- 
try w = — £ of the reduced equation. Finally, when expressed in terms of the 
rectifying coordinates v — X and Q = — ln(£) + X of w, Eq. (|4.6p reduces to 

-2c 3 v(l + v) -c 3 v 2 Q' = 0, (4.7a) 

with easy solution 

Q = A + — +v\n(v), (4.7b) 

c 3 v 

where A is an arbitrary constant of integration. However, trying to revert back to 
the solution w(z) to the original equation l|3.22a[) only leads to an implicit expression 



I + ProductLog 



exp(-l-^)) 



(4.4) 
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for w(z) of the form 



C3Z ProductLog(/(z, w)) 
where 



= B (4.8a) 



f(z,w) = -2exp(-A/c 3 )(2 + c 3 zw)/(c 3 z), (4.8b) 

and where A and B are the two arbitrary constants of integration. Consequently, 
it is not obvious how to determine the corresponding solution u to the ZK equa- 
tion and its properties from Eq. (|4.8[) . Similar difficulties are encountered with all 
Type B reduced equations, which turn out to have only zero- or one-dimensional 
symmetry algebras, and which by essence are more complicated then the Type A 
ones we have considered. 



4.1. Linearizability. Our analysis in this section revealed that in general all non- 
linear reduced ODEs were not readily integrable. Most often this is due to the fact 
that their solutions is not expressible in terms of algebraic or elementary functions. 
Before attempting to solve nonlinear ODEs obtained by similarity reductions, it is 
customary to perform a painleve test [19l [20] that determines whether these equa- 
tions satisfy certain necessary conditions for having the so-called Painleve property, 
because Painleve-type equations are generally perceived to be easier to solve than 
general nonlinear ODEs. For a given differential equation, the Painleve property 
is the property that all solutions of the equation are free from movable essential 
singularities, and for second order ODEs a list of fifty canonical forms of all such 
equations exists [2Tj . 

However, since we have attempted to solve all the nonlinear reduced ODEs with- 
out performing any Painleve test on them, we will simply contempt ourselves with 
checking whether the nonlinear reduced equations obtained, each of which may be 
put into the form 

A = w" + F(z, w, w') = 0, (4.9) 

for a certain function F, are linearizable by a transformation of the form 

z = f(Q,r), w = g{Q,r). (4.10) 

This will give a better insight into the integrability of the ZK equation, because 
linear equations are by essence integrable. Indeed, in the modern era of integra- 
bility, the Painleve approach has had an immense success by leading in particular 
to the discovery of a wealth of new integrable systems, to the extent that various 
conjectures almost intimately associated the Painleve property with integrability 
[23] . However, it has become clear that although many linearizable equations (both 
ordinary and partial) possess the Painleve property, many other linearizable equa- 
tions do not possess this property [22l[23]. By a result of Lie, the linearization 
test is easy to perform. Indeed, Lie [24 showed that a necessary condition for an 
equation of the form (|4.9j) to be linearizable is that it be of the form 



A = w" + A(z, w)w' 3 + B(z, w)w' 2 + C(z, w)w' + D(z, w) = 0, (4.11) 



9 



that is, a polynomial in w' of degree at most 3. If for a given equation of the above 
form (|4.11j) we consider the two expressions 

*i(A) = 3A Z ^ Z ~ 2B Z , W + C WtW - 3(CA) Z + 3(DA) W + (B 2 ) z + 3AD W - B{C) W 

*a(A) = 3D WiW - 2C Z , W + B ZiZ - 3{DA) Z + 3{DB) W - (C 2 ) w - 3DA Z + CB Z , 

then Lie also showed [2~H [2"5] that an equation of the form (|4.11[) is linearizable if 
and only if 

*l(A)=0, and * 2 (A) = 0, (4.12) 
and he also provided a method for finding the linearizing transformations of the 
form dUTDJ). 

All the nonlinear ODEs obtained for the ZK equation by similarity reductions 
are clearly of the form (14. lip required for the necessary condition of linearization. 
However, we have found that none of them is linearizable in general, except in some 
few cases where linearization is possible for a special value of the set of parameters. 
Indeed, if for every equation A of the form (143)) we set *(A) = (*i(A), * 2 (A)) , 
then it is readily seen that in the case of Eq. (|3 . 1 1 [) for instance, we have ^(A) = 
if and only if ko = 0, and for this value of ko the equation reduces to 

12 i II o, 

W + WW = (J, 

which transforms into the linear equation Z" = for Z = w 2 . A similar obser- 
vation is made about other Type A equations containing parameters. Indeed, the 
test reveals that Eqs. ()3.14a|) (|3.14bj) . and (I3.22a|) are linearizable if and only if 
ci = 0, C2 = and C3 = 0, respectively, in which case the equation is already 
linear. However, in the latter case (|3.22ap . the equation is in reality not linearizable 
because the value C3 = is not allowed. 

On the other hand none of the Type A equations without parameters is lineariz- 
able. For instance for Eq. (|3.2p we have 

54z 72(-2 + wzY 



*(A) - , 

yiA + wz) 3 ' z(A + wz) 3 

Similarly, none of the Type B equations, generally represented by (|3.4p . (]3 . 13[) 
and (|3.2ip is linearizable. For instance, for Eq. (|3.13p it is found that 'I'i(A) = if 
and only if ko = 1/9, and for this value of fco we have 

^- / a \ 18z 4/5 

WA) = ; . 

V ; 25(w + z 14 /5)2 

For Eq. ^^^^A) = if and only if fc = -1/15, while * 2 (A) = if and only if 
fco = 1/21. The linearization test thus shows that all nonlinear ODEs obtained by 
similarity reductions of the ZK equation are essentially non linearizable. 

5. Concluding remarks 

Based on the classification list of all two-dimensional subalgebras of the (2 + 1)- 
dimensional ZK equation recently proposed in [7] , we have obtained for the first time 
in this paper all possible similarity reductions to ODEs of the ZK equation. In this 
way, new exact solutions of this equation which are larger in number and different 
from those obtained in [6, 3J have been obtained. In general, all nonlinear reduced 
equations were not readily integrable, and instead of performing the usual Painleve 
test on these equations, we have rather shown the nonlinear reduced equations are 
essentially non linearizable. 
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